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Abstract 

•A  nonlinear  scene  analysis  algorithm  is  studied  for  complex 
scenes  containing  realistic  targets  and  background  clutter.  Infrared 
and  visible  spectrum  light  images  are  processed.  The  algorithm  com¬ 
bines  the  Fourier  transform  phase  array  of  a  scene  with  the  Fourier 
transform  magnitude  of  a  template  to  create  a  new  image.  Clutter 
reduction  ability  and  target  recognition  capability  are  examined  in 
detail. 
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SCENE  ANALYSIS  -  APPLICATION  OF 


TWO  DIMENSIONAL  NONLINEAR  FILTERING 
FOR  TARGET  ENHANCEMENT  AND  RECOGNITION 

I  Introduction 

Machine  detection  and  recognition  of  a  target  in  a  complex  or 
cluttered  scene  is  a  major  concern  in  advanced  military  technology  (Ref 
1,2).  Areas  of  application  are  analysis  of  reconnaissance  photos,  wea¬ 
pon  delivery  display  systems  and  "smart"  munitions. 

Background 

A  target  acquisition  process  can  be  viewed  as  consisting  of 
three  main  procedures.  The  first  procedure  is  the  development  of  a 
representation  of  a  scene  that  can  be  used  as  input  data.  The  second 
procedure  takes  the  input  data  and  extracts  characteristic  features  and 
attributes.  The  final  procedure  is  the  classification  and  identifica¬ 
tion  of  the  extracted  features  (Ref  3).  This  paper  assumes  that  the 
first  procedure  is  well  defined  and  that  the  input  data  used  are  images 
formed  from  light  in  the  visible  spectrum,  or  images  formed  by  energy 
from  the  infrared  spectrum. 

The  major  problem  of  detection  and  recognition  of  targets  in 
cluttered  scenes  is  extracting  the  characteristic  features  needed  to 
classify  the  target.  The  target  can  be  embedded  in  a  virtually  infi¬ 
nite  variety  of  backgrounds  with  many  of  these  backgrounds  having  the 
same  general  features  and  characteristics  as  the  target.  This  causes 
the  final  procedure  of  target  classification  and  identification  to  be 
extremely  difficult.  No  general  solution  to  this  problem  is  known. 
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In  1980,  Moshe  Horev  developed  a  picture  correlation  algorithm 
to  detect  and  identify  an  object  embedded  in  a  cluttered  scene.  This 
detection  could  be  accomplished  without  prior  knowledge  of  the  objects’ 
size,  orientation  or  location  (Ref  4).  This  algorithm  combined  the 
feature  extraction  process  and  the  identification  process  in  a  template 
matching  scheme  that  is  accomplished  in  two  stages: 

1.  The  first  stage  begins  by  a  series  of  transfor¬ 
mations  that  places  the  template  and  scene  into  a 
special  "correlation"  plane.  A  correlation  is  per¬ 
formed.  The  location  of  the  peak  value  of  the 
correlation  provides  the  angle  of  rotation  and  the 
difference  in  scale  between  the  target  and  the 
template.  The  template  is  then  modified  so  that  it 
is  at  the  same  scale  and  angle  of  rotation  as  the 
target. 

2.  In  the  second  stage  the  phase  angle  array  from  the 
Fourier  transform  of  the  image  is  combined  with  the 
magnitude  array  of  the  Fourier  transform  of  the 
modified  template.  Throughout  the  remainder  of  this 
thesis  this  process  will  be  referred  to  as  the  PIMT 
(Phase  of  Image,  Magnitude  of  Template)  process. 

An  inverse  Fourier  transform  of  the  combined  array 
is  computed.  The  result  is  a  modified  image,  in 
which  the  target  is  enhanced  while  the  clutter  is 
suppressed.  . 

Horev  used  the  algorithm  to  successfully  identify  targets  in 
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two  photographs  with  each  photograph  having  a  different  background  clut¬ 
ter.  The  significance  of  these  successes  is  the  apparent  ability  of  the 
algorithm  to  suppress  background  clutter.  In  a  separate  study  (Ref  5), 
infrared  pictures  were  processed  using  the  algorithm.  The  results  were 
unsuccessful.  No  detailed  study  has  been  made  of  the  reasons  of  the 
success  or  failure  of  the  algorithm. 

Problem 

The  purpose  of  this  thesis  is  to  investigate  the  effectiveness 
of  the  second  stage  of  the  Horev  algorithm  in  suppressing  background 
clutter. 

Equipment ,  Software  and  Data  Base 

The  image  processing  for  this  project  was  accomplished  in  the 
digital  signal  processing  laboratory  at  the  Air  Force  Institute  of 
Technology. 

Visible  spectrum  images  were  developed  in  the  laboratory  using 
a  system  consisting  of  a  CVC-1  Vidicon  camera,  video  moniter  and  a  NOVA 
2/10  minicomputer  interfaced  with  a  Cromemco  microcomputer.  The  images 
produced  were  256  X  256  pixel  arrays,  with  each  pixel  grey  scale  value 
ranging  from  0  to  15. 

The  infrared  data  base  consists  of  945  digitized  infrared 
images  provided  by  the  Air  Force  Armament  Laboratory,  Elgin  AFB,  Florida. 
The  images  are  in  TABILS  format  on  magnetic  tape.  Each  image  is  a 
100  X  100  pixel  array,  with  each  pixel  grey  scale  value  ranging  from 
0  to  1023.  Since  this  thesis  is  concerned  only  with  the  processing 
of  the  final  digitized  image,  a  detailed  description  of  the  data  base 
is  not  provided.  However,  a  detailed  description  of  the  TABILS  format 
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data  base  can  be  obtained  in  chapter  II  of  reference  5. 


All  computations  required  for  image  processing  were  accom¬ 
plished  on  a  Data  General  Eclipse  S/250  that  shares  disk  memory  with 
the  NOVA/Cromemco  system. 

The  photographs  contained  in  this  thesis  were  obtained  from  the 
NOVA/Cromemco  system  video  moniter.  A  complete  description  of  the  AFIT 
image  processing  facility  hardware  and  software  is  contained  in  refe¬ 
rence  6.  Software  unique  to  this  thesis  is  contained  in  the  appendix. 

Scope 

The  PIMT  process  was  used  on  filtered  and  unfiltered  images 
composed  of  different  scenes  and  targets.  A  threshold  operation  was 
performed  on  the  original  image  and  the  PIMT  image.  If  the  clutter  was 
reduced  and  the  target  enhanced  in  the  PIMT  image,  when  compared  to  the 
original  image,  the  process  was  considered  a  success.  This  study  does 
not  attempt  to  develop  a  rigorous  mathematical  analysis  of  the  non¬ 
linear  PIMT  process.  However,  a  theory  is  developed  and  presented 
based  on  a  study  of  the  results  of  the  process. 

Sequence  of  Presentation 

This  study  begins  with  a  discussion  of  the  PIMT  process 
theory  and  implementation.  Chapter  III  discusses  the  processing 
used  to  create  the  scene  and  template  images  used.  Chapter  IV  veri¬ 
fies  that  the  process  of  Chapter  II  (PIMT  Process)  obtains  the  same 
results  as  those  obtained  by  the  second  stage  of  the  Horev  algorithm. 
Chapters  V  and  VI  examine  the  effect  of  filtering  on  the  PIMT  process. 
Chapter  VII  examines  the  effect  of  size  and  rotation  variations  and 
Chapter  VIII  examines  the  discrimination  ability  of  the  process. 
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II  PIMT  Process 


Theory 

The  PIMT  process  is  the  combination  of  the  phase  array  from  the 
Fourier  transform  of  the  scene  combined  with  the  magnitude  array  of  the 
Fourier  transform  of  the  template.  Let  the  intensity  functions  of  the 
template  image  and  scene  image  be  designated  as  t(x,y)  and  s(x,y),  respec¬ 
tively.  Let  the  2-D  I>FT  be  designated  as  in  equations  1  and  2. 


F{ t(x,y) }  =  T(c,n)  (1) 

F{s(x,y) }  =  S(?,n)  (2) 

The  Fourier  transform  is  a  complex  function  and  can  be  represented  as: 

TU.n)  =  R  (5,n)  +  j  l  a>n)  (3) 

s(s,n)  =  R  (c,n)  *  j  I  (?»n)  (4) 


where  R(?,n)  and  I(C,n)  are  the  real  and  imaginary  components  of  the 
transformation.  Using  Euler's  equation,  T(t,n)  and  S(C,n)  can  be  repre¬ 
sented  in  terms  of  their  magnitude  and  phase  spectrum  as  shown  below. 
T(s,n)  *  |T(s,n)  |exp[ j<j>t(S»u)  ]  (5) 

S(c,n)  =  js(c,n)  |exp[j<t>s(*;,n)  ]  (6) 

It  should  be  noted  that  T(c,n)  and  S (^ , n )  are  each  represented  as  a  dis¬ 
crete  256  X  256  complex  array.  Every  point  in  that  array  has  a  magnitude 
and  phase  associated  with  that  point.  The  result  of  the  PIMT  process  can 


be  designated  as  shown  in  equation  7. 

P(C,n)  =  |T(c,n)  |exp[j4>gU,n)]  (7) 

Looking  at  a  particular  point  C=a,  n=b  then 

T(a,b)  =  |T(a,b)  |exp[j<t>t(a,b)]  (8) 

S(a,b)  =  |S(a,b) | exp [ (a,b) ]  (9) 

P(a,b)  =  |T(a,b)  |exp[j4>s(a,b)]  (10) 
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The  point  (a,b)  in  the  Fourier  (spatial  frequency)  domain  repre¬ 
sents  a  location  on  the  spatial  frequency  plane.  If  the  scene  consists 
of  background  clutter  plus  target,  then  the  phase  component  of  the  fre¬ 
quency  spectrum  is 

<t>g(a,b)  =  4>t  (a,b)  +  0(a,b)  (11) 

where  9(a,b)  is  the  phase  change  due  to  the  background  clutter  or  noise 
present  in  the  scene.  If  the  noise  component,  6(a,b)  is  zero  or  very 
small  for  all  (a,b),  then 

P(a,b)  *  T(a,b)  (12) 

and  the  image  of  the  target  would  be  essentially  unaffected  in  the  PIMT 
image.  If,  however,  the  term  0(a,b)  7s  0  for  all  (a,b),  then  the  target 
information  will  be  distorted  in  the  PIMT  image.  Similarly,  if  the  tem¬ 
plate  magnitude,  |T(a,b) j ,  is  less  than  the  scene  magnitude,  js(a,b)|, 
for  any  (a,b),  then  the  target  energy/clutter  energy  ratio,  or  the  signal 
to  noise  ratio  should  be  increased  causing  a  reduction  in  the  background 
clutter.  Thus,  the  PIMT  process  can  be  seen  to  be  optimized  if  the  spa¬ 
tial  frequencies  of  the  clutter  differ  significantly  from  the  frequencies 
of  the  target  making  the  PIMT  process  both  scene  and  target  dependent. 

This  will  be  demonstrated  in  the  following  sections. 

Implementation 

The  PIMT  process  was  implemented  as  shown  in  the  flow  chart  of 

Figure  1. 

The  first  step  involved  computing  two  dimensional  fast  Fourier 
transforms  (2-D  DFT)  for  both  the  template  and  scene  images.  If  filtering 
was  used,  it  occurred  in  the  spatial  frequency  domain  as  shown  in  the 
processing  steps  designated  A  and  B  of  Figure  1.  The  PIMT  process  was 
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then  performed  to  combine  the  Fourier  transform  phase  array  of  the  image 
with  the  Fourier  transform  magnitude  array  of  the  template.  An  inverse 
2-D  DFT  was  then  performed  on  the  result  of  the  PIMT  process  to  form  the 
PIMT  image.  For  the  remainder  of  this  thesis,  the  forward  and  inverse 
2-D  DFT’s  will  be  considered  an  integral  part  of  what  will  be  referred  to 
as  the  PIMT  process. 
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Ill  Image  Processing 


The  image  processing  used,  in  addition  to  the  PIMT  process,  were 
the  processes  used  to  create  the  images  and  those  used  to  filter  the 
created  images.  This  chapter  will  contain  a  discussion  of  the  procedures 
used  to  create  the  scene  and  template  images.  A  discussion  of  the  filter¬ 
ing  processes  will  be  presented  in  Chapter  VI. 

Infrared  Spectrum 

All  processing  of  the  infrared  images  (including  filtering) 
utilized  the  full  1024  grey-scale  levels  provided  by  the  data  base.  The 
photographs  shown  in  this  section  and  the  remainder  of  this  thesis  were 
obtained  from  the  NOVA/Cromemco  system  video  monitor.  Since  this  system 
only  uses  a  16  level  grey  scale,  a  linear  scaling  transform  (Ref  7:161) 
was  performed  as  a  final  step  to  allow  display  of  the  infrared  images. 

A  template  of  a  tank  was  created  from  the  first  IR  scene  of  the 
data  base.  This  scene  is  shown  in  Figure  2. 

The  procedure  used  to  create  the  template  was  an  edge  search 
operation  as  shown  in  Figure  3.  The  operation  starts  at  the  right  edge 
of  the  scene  (pixel  Aq)  and  computes  the  difference  between  two  adjacent 
pixels  in  the  row  (Aj^  -  AN_^) .  If  the  difference  is  below  or  equal  to 
a  certain  threshold  (T),  the  difference  is  then  calculated  and  tested 
for  the  next  two  pixels.  This  continues  until  the  difference  exceeds 
the  threshold  value.  When  this  occurs,  all  of  the  pixels  in  the  row  from 
Aq  to  A^  are  set  equal  to  zero  completing  the  operation.  The  same  pro¬ 
cedure  was  applied  to  the  top,  left,  and  bottom  edges.  The  resulting 
template  is  shown  in  Figure  4.  This  procedure  worked  well  in  the  scene 


Figure  4.  Infrared  Tank  Template 


of  Figure  2  because  there  were  no  high  energy  spots  (hotspots)  in  the 
background  and  the  tank  was  significantly  "hotter"  than  the  rest  of  the 
scene.  The  procedure  is  not  a  general  one  and  it  is  likely  that  it  would 
fail  when  applied  to  most  scenes. 

A  three  by  three  image  array  was  created  as  a  master  scene  to 
test  the  PIMT  process.  The  nine  infrared  images,  including  the  scene  of 
Figure  2,  were  normalized  to  insure  that  no  particular  image  would  be 
"hotter"  than  any  other  particular  image  used  in  the  master  scene.  The 
nine  images  were  combined  to  produce  the  master  scene  shown  in  Figure  5. 

Visible  Spectrum 

Visible  spectrum  images  were  created  in  the  laboratory  from 
black  and  white  photographs.  Each  photograph  was  digitized  into  a  256  X 
256  pixel  array.  Random  noise  introduced  during  digitization  was  reduced 


Figure  5.  Master  Scene  of  Infrared  Images 


using  frame-to-frame  averaging  (Ref  7:206).  The  averaged  image  was  then 
tested  for  the  minimum  and  maximum  grey  scale  level.  A  linear  scaling 
transform  (Ref  7:161)  was  used,  when  necessary,  to  fully  utilize  the  16 
grey  levels  available  with  the  NOVA/Cromemco  system. 

Templates  were  created  by  first  segmenting  the  desired  object 
into  a  rectangular  window  surrounded  by  a  zero  grey-scale  level  back¬ 
ground.  A  "trial-by-error"  thresholding  was  then  done  on  the  window 
until  only  the  desired  template  remained.  Figure  6  shows  a  black  and 
white  photograph  prior  to  digitization.  Figure  7  shows  the  result  of  the 
digitization  and  noise  reduction  process  and  Figure  8  is  a  completed 


template. 


I 

Figure  8.  Final  Template 

The  template  could  be'  moved  to  any  location  within  the  array  and 
combined  with  any  background  scene.  This  process  guaranteed  that  the 
target  in  the  scene  would  be  the  same  scale  and  rotation  as  the  template 
used,  in  accordance  with  the  original  assumption. 
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IV  Verification  of  the  PIMT  Process 


The  stated  purpose  of  this  thesis  is  to  investigate  the  effec¬ 
tiveness  of  the  second  stage  of  the  Horev  algorithm  in  suppressing  back¬ 
ground  clutter.  Two  samples  of  results  are  presented  in  this  chapter  to 
demonstrate  that  background  clutter  reduction  produced  by  the  PIMT  process 
is  similar  to  that  obtained  using  the  Horev  algorithm.  Since  the  scenes 
and  templates  used  by  Horev  were  not  available  for  this  test,  success  is 
based  solely  on  the  similarity  of  the  clutter  reduction  obtained. 

Procedure 

PIMT  processing  was  the  only  procedure  used  to  obtain  the  results 
contained  in  this  chapter.  The  scenes  used  are  shown  in  Figures  9A  and 
9B.  The  template  used  for  both  scenes  is  shown  in  Figure  10. 


Figure  9A.  First  Scene  Image 
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Template 


Results 


The  final  PIMT  images  are  shown  in  Figures  11A  and  11B.  It  can 
be  seen  by  comparing  Figures  9A  and  11A  and  Figures  9B  and  11B  that  a 
significant  reduction  in  background  clutter  occurs  as  a  result  of  the  PIMT 
process.  This  does  demonstrate  that  the  PIMT  process  reduces  clutter 
similar  to  the  clutter  reduction  obtained  by  the  Horev  algorithm.  Further 
analysis  of  the  images  shown  in  Figures  11A  and  11B  will  be  discussed  in 
chapter  IX. 


I 

Figure  11A.  PIMT  Image  of  the  Scene  of  Figure  9A 
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I  Figure  11B.  PIMT  Image  of  the  Scene  of  Figure  9B 
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V  Small  Window  Processing 


Horev  hypothesized  that  the  PIMT  process  works  best  when  a 
small  window  is  used  (Ref  4:87).  The  theory  is  that  there  will  be  a  low 
clutter-energy  to  target-energy  ratio,  since  most  of  the  scene  will  con¬ 
sist  of  the  target.  The  infrared  images  were  already  in  a  small  window 
format  and,  therefore,  used  to  test  Horev' s  hypothesis. 

Procedure 

The  PIMT  images  were  created  using  the  template  shown  in  Figure 
4  and  scenes  shown  in  Figures  2,  13A  and  14A.  Correlations  were  per¬ 
formed  between  the  template  image  and  both  the  PIMT  image  and  the  original 
scene  image.  The  two  correlations  were  then  compared  to  see  if  the  PIMT 
processing  resulted  in  a  higher  correlation  between  the  processed  scene 
and  the  template. 

Results 

Figure  12  is  the  PIMT  image  produced  using  the  template  of 
Figure  4  and  the  scene  of  Figure  2.  Figures  13A  and  14A  are  original 
scene  images  and  Figures  13B  and  14B  are  the  resulting  PIMT  images  re¬ 
spectively. 

It  should  be  noted  that  even  with  a  small  window  there  is  a 
tendency  for  the  PIMT  image  to  be  merely  a  distorted  reproduction  of  the 
original  scene.  A  comparison  of  the  correlations  between  the  PIMT 
images  and  the  template,  and  of  the  correlations  between  the  original 
scenes  and  the  template  shows  that  the  correlations  using  the  original 
scenes  worked  better  in  target  detection.  In  either  case,  there  was  no 
significant  difference  between  the  two  operations. 
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Figure  13B.  PIMT  Image  of  Scene  in  Figure  13A 


Figure  14A.  Scene  Containing  No  Target 


Figure  14B.  PIMT  Image  of  Scene  in  Figure  14A 
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VI  Filtering 


Various  filters  were  used  in  conjunction  with  the  PIMT  process. 
All  filtering  was  accomplished  in  the  spatial  frequency  domain  as  shown 
in  the  processing  steps  designated  as  A  and  B  in  Figure  1.  This  chapter 
will  briefly  describe  the  filters  used  on  the  scene  images  and  discuss 
the  results  on  the  PIMT  image. 


Filters 


The  filters  used,  as  described  in  this  section,  are  a  function 
of  a  radial  frequency  f.  This  radial  frequency  describes  a  radius  in  the 
frequency  domain  starting  from  the  centered  DC  component  to  f.  The  follow¬ 
ing  filters  were  used: 


H(f)  =  1  (no  filtering) 
H(f)  «  f 
H(f )  =  1/f 
H(f )  =  f2 


H(f) 


(f  S  12) 


H(f )  = 


6 


(fx  1  f  1  f2) 
(0  elsewhere) 


(13) 

(14) 

(15) 

(16) 

(17) 

(18) 


The  filters  described  in  Equations  14  and  16  were  used  to 
duplicate  the  effect  of  the  transformations  used  by  Horev  in  his  first 


stage  of  processing  (Ref  4:40).  The  filter  of  equation  15  is  essentially 


a  low-pass  filter  and  was  used  to  investigate  the  effects  of  low-pass 
filtering  on  the  PIMT  process.  The  band-pass  filter  shown  in  equation 
18  was  used  to  accomplish  a  frequency  analysis  of  the  infrared  image  of 
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Figure  5.  It  was  used  to  find  a  frequency  range  in  which  the  PIMT  re¬ 
sults  would  be  more  successful.  Photographs  of  these  filtered  scenes  and 
resulting  PIMT  image  are  included  at  the  end  of  this  chapter. 

Three  other  filter  type  operations  were  used  in  an  attempt  to 
improve  the  PIMT  process.  These  operations  compared  the  scene  with  tem¬ 
plate  to  try  to  emphasize  those  frequency  components  of  the  template 
which  were  much  stronger  than  corresponding  frequency  components  of  the 
image.  Anew  template,  T'(5,n)  was  created  as  shown  in  Equations  19 


through  21  below. 


(|TU,n)|-lsu,n)|  if  |T(c,n)|-|s(c,n)|>0 
|t' (c,n)  |  =  (19) 

.  0  Elsewhere 


|t'U,ti)|  = 


|t(s, n) | 

|s(c,n)| 


|T'(c,n)| 


|T(c,n)l 

|T(c,n)|-|s(c,n)| 


Equation  21  was  used  because  Horev  predicted  that  it  would  automatically 
yield  "the  'best*  set  of  features"  (Ref  4:102). 


Results 


Figures  15  through  47  show  the  results  of  the  various  filters 


on  infrared  images.  It  can  be  seen  that  none  of  the  filters  caused  sig¬ 
nificant  improvement  in  the  PIMT  process.  Similar  results  were  obtained 
using  isnag  from  the  visible  spectrum. 

These  results  are  best  interpreted  using  Equation  7.  In  the 
analysis  in  Chapter  II,  it  was  proposed  that  the  clutter  reduction 
largely  came  about  when,  for  any  particular  4»a  and  ri=b,  | T (a , b) |  was 
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very  small  or  equal  to  zero.  This  causes  the  corresponding  phase  com¬ 
ponent  4>  (a,b)  to  be  ineffective.  This  would  reduce  the  non-template 
s 

frequency  components  of  the  scene  image.  While  the  corresponding  tem¬ 
plate  frequency  components,  i.e.  the  target,  would  be  unaffected.  Since 
the  filtering  caused  no  improvement  in  the  PIMT  process  there  are  at 
least  three  possible  conclusions.  One,  the  right  filter  function  had 
not  been  applied.  Two,  the  scene  and  the  template  have  similar  frequency 
magnitude  spectrums.  Three,  the  PIMT  process  is  independent  of  the  tem¬ 
plate  frequency  magnitude  spectrum. 

Reference  8  and  Chapter  VIII  of  this  thesis  present  experimen¬ 
tal  evidence  that  for  target  identification,  the  PIMT  process  is  indepen¬ 
dent  of  the  template  magnitude.  Although  no  detailed  comparison  was 
made,  examination  of  the  2-D  DFTs  of  the  scenes  and  the  templates  have 
shown  that  both  have  similar  spectrums.  Consistent  failure  of  spectral 
filtering  techniques  when  applied  to  complex  scenes  was  the  basis  for 
assuming  similar  spectrums  at  the  start  of  this  research.  Given  the 
evidence  supporting  the  last  two  conclusions  it  is  the  authors’  opinion 
that  there  is  no  "right"  filter. 
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Figure  17.  PIMT  Image  of  Figure  16 


Figure  18.  Filtered  Masterscene:  H(f)  *  1  for  1.5  f  525,  0  Elsewhere 
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Figure  19.  PIMT  Image  of  Figure  18 
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Figure  20.  Filtered  Masterscene: 
H(f)  ■  1  for  10 <  f  <25,  0  Elsewhere 
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Figure  22.  Filtered  Masterscene: 
[(f)  “  1  for  10  if <  50,  0  Elsewhere 


Figure  24.  Filtered  Masterscene: 
H(f)  *  1  for  25<fl50,  0  Elsewhere 
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Figure  27.  P1MT  Image  of  Figure  26 


Figure  28.  Filtered  Masterscene: 
H(f)  =  1  for  100  <  f  <  256,  0  Elsewhere 
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Figure  37.  Masterscene  with  Different  Target  (Center) 


37 


$ 


Figure  41.  Filtered  Scene  of  Figure  40:  H(f)  =  f 


Figure  43.  Scene  with  Different  Target  (2) 


U  I 

Figure  44.  PIMT  Image  of  Figure  43:  H(f)  =  1 
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Figure  45.  PIMT  Image  of  Figure  2:  H(f)  =  1/f 


Figure  46.  PIMT  Image  of  Figure  2:  H(f)  =  f 


VII  Si2e  and  Rotation  Studies 


The  next  two  chapters  present  the  results  of  studies  made  to 
determine  the  discrimination  abilities  of  the  PIMT  process.  These  studies 
are  concerned  with  how  much  size  and  rotational  deviation  of  the  target 
from  the  template  being  used  is  allowed  by  the  PIMT  process.  This  chap¬ 
ter  presents  the  results  of  target  size  and  rotation  deviation  studies. 

The  next  chapter  presents  the  results  on  how  well  the  PIMT  process  can 
discriminate  one  target  from  another  (i.e.  a  truck  from  a  tank). 

Procedure 

The  only  procedure  used  for  these  studies  was  the  PIMT  process. 
No  filtering  of  the  scene  or  template  images  was  accomplished.  The  scene 
and  template  images  were  created  from  the  black  and  white  photograph 
shown  in  Figure  6.  The  photograph  was  digitized  by  the  camera  located 
at  various  distances  from  the  photograph,  resulting  in  images  of  dif¬ 
ferent  size.  Rotated  images  were  created  by  rotating  the  photograph  in 
five  degree  increments  from  zero  to  ninety  degrees,  keeping  the  size  of 
the  photograph  constant. 

Figure  48  shows  the  template  image  used  for  the  first  test. 
Figures  49  and  50  show  the  scene  image  and  the  resulting  PIMT  image. 

As  can  be  seen,  most  of  the  details  of  the  original  scene  are  retained 
in  the  PIMT  image  even  though  the  template  used  was  approximately  1/10 
smaller. 

For  the  next  test,  the  template  of  Figure  51  and  the  scene  of 
Figure  53  below  were  processed.  The  target  in  this  scene  is  1/3  larger 
and  rotated  45  degrees  from  the  template  used.  For  comparison,  a  PIMT 
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Figure  48.  Template  Used  for  Size  Study 
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Figure  49.  Scene  Image  Used  for  Size  Study 
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Figure  50.  PIHT  Image  Used  for  Size  Study 


image  was  also  created  using  the  same  scene  (Figure  53)  and  the  tem¬ 
plate  of  the  rotated  target  in  the  scene  as  shown  in  Figure  52.  The 
resulting  PIMT  images  are  shown  in  Figures  54  and  55,  respectively. 

Again,  it  can  be  seen  in  Figure  54  that  the  clutter  is  redu¬ 
ced  and  the  target  is  visible  using  a  target  three  times  larger  than 
the  template  used. 

Discussion 

This  study  demonstrates  that  although  optimal  PIMT  clutter  re¬ 
duction  and  potential  target  identification  occurs  when  the  template  is 
the  same  size  and  rotation  as  the  target  in  the  scene,  changes  in  size 
and  rotation  still  result  in  successful  clutter  reduction.  The  target 
is  distorted,  but  not  to  the  point  that  potential  target  identification 
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would  not  occur.  It  can  be  concluded  that  severe  size  and  rotation 
deviations  are  allowed  by  the  PIMT  process. 
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Figure  51.  Template  of  Reduced  and  Unrotated  Target 
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Figure  53.  Scene  Image  with  Rotated  Target 
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Figure  54.  PIMT  Image  of  Unrotated  Template 


Figure  55.  PIMT  Image  of  Rotated  Template 

I 
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VIII  Discrimination  Studies 


This  chapter  is  a  report  on  a  study  made  to  determine  the  dis¬ 
crimination  ability  of  the  PIMT  process.  The  intent  was  to  see  if  the 
process  could  differentiate  between  a  truck  and  a  tank  or  any  other  wea¬ 
pon  system. 

Procedure 

The  PIMT  process  was  used  with  no  filtering.  The  scene  image 
used  for  processing  is  shown  in  Figure  56.  Templates  were  created  of 
each  of  the  objects  and  PIMT  images  were  obtained.  Figure  57  shows  the 
template  created  from  the  top,  middle  tank  and  figure  58  is  the  resulting 
image . 


I 

Figure  56.  Scene  for  Discrimination  Test 
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For  all  templates,  the  PIMT- images  showed  an  equal  energy 
distribution  among  the  objects  similar  to  that  shown  in  the  photograph 
above  (i.e.  the  same  PIMT  image  no'matter  which  -template  was  used). 

It  was  suspected  that  a  similar  result  could  be  obtained  from 

\ 

a  "dummy  template",  a  template  containing  no  details  such  as  gun  barrels 
or  turrets  (Ref  8).  A  template  was  created  consisting  of  a  rectangle  of 
constant  grey  scale  levej.  that  had  the  approximate  dimensions  of  the 
objects  in  the  scene'image.  The  resulting  PIMT  image  is  shown  in  the 

v  N“ 

photograph  below. 


Figure  59.  PIMT  Image  Using  Dummy  Template 


Results 

The  results  of  this  study  show  that  the  PIMT  process  has  little 
or  no  discrimination  ability. 
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IX  Conclusion 


The  purpose  of  this  thesis  was  to  investigate  the  second  stage 
of  the  Horev  algorithm  in  suppressing  background  clutter  in  infrared  and 
visible  spectrum  images.  This  thesis  presents  representative  results 
obtained  from  processing  approximately  100  infrared  spectrum  images  and 
80  visible  spectrum  images. 

Conclusion 

It  can  be  concluded  that  the  second  stage  of  the  Horev  algo¬ 
rithm,  the  PIMT  process,  does  suppress  background  clutter.  This  was 
verified  by  the  results  presented  in  Chapter  IV. 

A  theory  is  presented  in  Chapter  II  that  demonstrates  the  PIMT 
process  suppresses  background  clutter  by  eliminating  or  reducing  the  fre¬ 
quency  components  contained  in  the  scene  image  which  are  not  contained 
in  the  template  image.  The  results  presented  in  Chapters  IV  through  VI 
demonstrate  that  the  clutter  reduction  is  highly  dependent  on  the  scene 
and  template  images.  This  dependence  is  best  demonstrated  by  the  two 
"successful"  images  presented  in  chapter  IV.  An  analysis  of  the  template 
shows  that  most  of  the  template  energy  will  be  concentrated  in  the  narrow 
band  of  spatial  frequencies  required  to  produce  the  horizontal  "lines"  in 
the  truck.  A  similar  analysis  of  the  scene  shows  that  the  energy  will  be 
spread  over  a  wide  range  of  frequencies.  It  can  be  seen  in  the  resulting 
PIMT  image  that  the  spatial  components  retained  from  the  original  scene 
image  are  those  components  of  the  narrow  band  of  frequencies  of  the  tem¬ 
plate.  Looking  back  again  at  the  original  scene,  it  can  be  seen  that 
most  of  the  energy  in  this  narrow  band  of  frequencies  is  concentrated  at 
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the  target.  This  combination  of  favorable  factors  results  in  the  suc¬ 
cessful  PIMT  image. 

Chapter  V  demonstrates  that  when  the  potential  target  is  the 
largest  component  of  the  scene,  the  PIMT  process  is  no  better  than  a 
straight  correlation  operation.  This  causes  doubt  as  to  the  usefullness 
of  the  PIMT  process  as  an  initial  stage  of  automatic  pattern  recognition. 
Chapter  VI  demonstrates  that  filtering  causes  no  improvement  of  the  scene 
image.  Chapter  VII  demonstrates  that  optimal  target  detection  and 
clutter  reduction  occurred  when  the  template  is  the  same  size  and  rota¬ 
tion  as  the  target  in  the  scene  image.  It  is  also  demonstrated  that 
severe  size  and  rotation  deviations  of  the  template  to  the  target  are 
allowed  by  the  PIMT  process.  Finally,  it  is  shown  in  Chapter  VII  that 
the  PIMT  process  has  little  or  no  discrimination  ability,  making  it  more 
of  a  "blob"  detection  process  than  a  pattern  recognition  process. 

The  results  obtained  in  this  study  are  consistent  to  those 
obtained  by  Horev  (Ref  4)  and  Fadem  and  Walters  (Ref  5).  It  is  concluded 
that  the  PIMT  process  is  unsuccessful  as  a  single  process  automatic 
pattern  recognition  algorithm  due  to  its  high  dependence  on  the  scene 
and  template  images  used. 
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Appendix 


Computer  Programs 

Tt.e  following.,  pages  are  the  computer  programs  used  that  are 
unique  to  this  thesis.  This  software  is  designed  to  be  functional  and 
is  not  optimized  for  speed  or  minimum  code. 
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PROGRAM  READTAPE 


THIS  PROGRAM  READS  DATA  FOR  ONE  PICTURE  FROM 
THE  CONVERTED  TABILS  FORMAT  TAPE,  AND  STORES  THE 
DATA  IN  A  FILE  CALLED  PICTURE. 


DIMENSION  10(256) 

INTEGER  COMMAND, B( 16) 

COMMAND=0 

CALL  INIT( wMT0  n , 0 , 1 ) 

IFCI. NE. 1 )TYPE  " INIT=  " , I 
CALL  DFILW( "PICTURE”, I) 

CALL  CFILW( "PICTURE", 3, 102, I) 

IF(I.NE.I)  TYPE  "CFILWs " , I 
CALL  OPEN ( 2 , "PICTURE" , 3 ,1,120) 

IF(I.NE.I)  TYPE  ”OPEN= " , I 
TYPE  "THE  TAPE  CONTAINS  6  FILES  WITH  40  PICTURES" 
TYPE  "PER  FILE.  ENTER  FILE  NUMBER , 1 , 2 , 3 , 4 , 5 , OR , 6 " 
ACCEPT  IFILE 

TYPE  "ENTER  PICTURE  NUMBER  1 ,2 ,3 , . . .37 *38, 39 , OR,40" 
ACCEPT  NUM 

POSITION  THE  TAPE  TO  THE  DESIRED  FILE 

GOTO( 10, 20, 30, 40, 50, 60) , IFILE 
10  CALL  MT0PDO  ,"MT0:0", 0,1) 

GOTO  70 

20  CALL  MT0PD(1 ,"MT0:1", 0,1) 

GOTO  70 

30  CALL  MT0PD(1 ,"MT0:2", 0,1) 

GOTO  70 

40  CALL  MT0PD(1 ,"MT0:3",0,I) 

GOTO  70 

50  CALL  MT0PD(1 ,"MT0:4", 0,1) 

GOTO  70 

60  CALL  MT0PD(1 ,"MT0:5", 0,1) 

70  CONTINUE 

IF(NUM.EQ.I)  GOTO  101 

SET  TAPE  TO  DESIRED  PICTURE.  CALL  TO  WORD  TELLS 
MTDIO  TO  SKIP  RECORDS.  NUM*102  IS  NUMBER  OF 
RECORDS  SKIPPED.  SEE  MTDIO  COMMAND  IN  FORTRAN  IV 
MANUAL. 

CALL  WORD( COMMAND, 0,0,1 ,1,0, 0,0, 0,0, 0,0, 0,0, 0,0,0) 

COMMAND= COMMAND+( NUM-1 ) «102 

CALL  MTDI0(1 , COMMAND , 10 , ISTAT, I , IRC) 

TYPE  "RECORDS  SKIPPED: " , IRC 
TYPE  " ISTAT: ", ISTAT 

IF(I.NE.I)  TYPE  "MTDI01:", I, "ISTAT:", ISTAT 
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101  COMMANDS 

DO  100  K= 1 ,102 
KK=K-1 

READ  A  RECORD  OFF  TAPE  INTO  ARRAY  10 

CALL  MTDI0( 1 , COMMAND , 10 , ISTAT , I , IRC) 

IF(I.NE.I)  TYPE  "MTDI0= " ,  I , "ISTAT= " , ISTAT 
TYPE  “RECORD", K, "  ",  IRC, "WORDS  READ" 

DELETE  SOME  GARBAGE  INSERTED  BY  MTDIO  COMMAND 

DO  88  JK=207,255 
88  IO(JK)=0 

INSERT  CARRIAGE  RETURNS  EVERY  128  CHARACTER, 

THE  CALL  TO  WORD28  SETS  THE  LAST  8  BITS  TO  00001101 
(BINARY  FOR  CR). 

DO  110  J=6M  ,256  ,64 
CALL  WORD28C IO( J) ,0 ,0 ,0 ,0 , 1  , 1  ,0 , 1 ) 

11U  CONTINUE 

WRITE  A  BLOCK  (ARRAY  10)  ONTO  FILE  NAMED  PICTURE 

CALL  WRBLK(2,KK,I0,1 ,IW) 

IF(IW.NE.I)  TYPE  "WRBLK= " , IW 
100  CONTINUE 

CALL  CLOSE( 1 , IC) 

IF(IC.NE.I)  TYPE  "CLOSE=",IC 
CALL  RLSE( "MTO " , IM) 

IF(IM.NE.I)  TYPE  "RLSE=",IM 
STOP  "ONE  PICTURE  READ" 

END 
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PROGRAM  WORD 


Subroutine  WORD  sets  the  bits  in  a  16  bit  word. 
If  B1 5=1  bit  15  is  set,  if  B15=0  bit  15  is  cleared. 


SUBROUTINE  WORD( I , B1 5 , B1 4 , B1 3 , B1  2 , B1 1  ,B10,B9,B8, 
+  B7,B6,B5,B4,B3,B2,B1 ,B0) 

INTEGER  BO,B1 ,B2,B3,B4,B5,B6,B7,B8,B9,B10,B11  , 

+  B12,B13,B14,815 

IF( BO . EQ. 0) GOTO . 1 0 
CALL  BSET( 1,0) 

GOTO  11 

10  CALL  BCLR(  1,0) 

11  IF(B1 .EQ.O)  GOTO  20 

CALL  BSET( 1,1) 

GOTO  21 

20  CALL  BCLR( 1,1) 

21  IF( B2 . EQ. 0 )  GOTO  30 

CALL  BSET( 1,2) 

GOTO  31 

30  CALL  BCLR( 1,2) 

31  IFCB3.EQ.0)  GOTO  40 

CALL  BSET( 1,3) 

GOTO  41 

40  CALL  BCLR(I,3) 

41  IF( B4 . EQ. 0 )  GOTO  50 

CALL  BSET(W0RD, 4) 

GOTO  51 

50  CALL  BCLR  (1,4) 

51  IF(B5.EQ.O)  GOTO  60 

CALL  BSET(I,5) 

GOTO  61 

60  CALL  BCLR( 1,5) 

61  IF(B6.EQ.0)  GOTO  70 

CALL  BSET( 1,6) 

GOTO  71 

70  CALL  BCLR( 1,6) 

71  IF( B7 .EQ.O)  GOTO  80 

CALL  BSET(I,7) 

GOTO  81 

80  CALL  BCLR( 1,7) 

81  IF(B8.EQ.0)  GOTO  90 

CALL  BSET( 1,8) 

GOTO  91 

90  CALL  BCLR( 1,8) 

91  IF( B9. EQ. 0)  GOTO  100 

CALL  BSET(I,9) 

GOTO  101 

100  CALL  BCLR(I.9) 

101  IF( BIO . EQ.O)  GOTO  110 
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CALL  BSET( 1,10) 

GOTO  in; 

110  CALL  BCLR(I.IO) 

m  IF( B1 1 .EQ.O)  GOTO  120 

CALL  BSET( 1,11) 

GOTO  121 

120  CALL  BCLR(I,11) 

121  IFCB12.EQ.0)  GOTO  130 

CALL  BSET( 1,12) 

GOTO  131 

130  CALL  BCLR( 1,12) 

131  IF(B13.EQ.O)  GOTO  140 

CALL  BSET(I,13) 

GOTO  141 

140  CALL  BCLR(I.13) 

141  IFCB14.EQ.0)  GOTO  150 

CALL  BSET( 1,14) 

GOTO  151 

150  CALL  BCLR( 1,14) 

151  IF( B1 5 . EQ. 0 )  GOTO  160 

CALL  BSET(I,15) 

GOTO  161 

160  CALL  BCLR( I , 1 5 ) 

161  CONTINUE 
RETURN 
END 


o  o  o  o  o  o  o 


PROGRAM  WORD28 


Subroutine  W0RD28  sets  the  last  8  bits  in  a  16  bit 
word.  If  B0=0  bit  0  will  be  cleared.  If  B0  =  1  bit  0 
will  be  set. 


SUBROUTINE  WORD28 ( IWORD , B7 , B6 , B5 , B4 , B3 , B2 , B1 ,B0) 
INTEGER  B0,B1 ,B2,B3,B4,B5,B6,B7,B8,B9,B10,B11 ,B12, 
+  B13,B14,B15 

IF( BO . EQ. 0)G0T0  10 
CALL  'BSET(  IWORD, 0) 

GOTO  11 

10  CALL  BCLR( IWORD, 0) 

11  IF(B1 . EQ.O)  GOTO  20 
CALL  BS£T( IWORD, 1 ) 

GOTO  21 

20  CALL  BCLR( IWORD, 1) 

21  IF( B2 . EQ. 0 )  GOTO  30 
CALL  BSET( IWORD, 2) 

GOTO  31 

30  CALL  BCLR( IWORD, 2) 

31  IF( B3 . EQ.O)  GOTO  40 
CALL  BSET( IWORD, 33 
GOTO  41 

40  CALL  BCLR( IWORD , 3 ) 

41  IF( B4 .EQ.O)  GOTO  50 
CALL  BSE7( IWORD, 4 ) 

GOTO  51 

50  CALL  BCLR  (IWORD, 4) 

51  IF(B5.EQ.O)  GOTO  60 
CALL  BSET( IWORD, 5 ) 

GOTO  61 

60  CALL  BCLR( IWORD, 5) 

61  IF(B6 .EQ.O)  GOTO  70 
CALL  BSET(IW0RD,6; 

GOTO  71 

70  CALL  BCLR( IWORD , 6 ) 

71  IF( B7 . EQ. 0)  GOTO  80 
CALL  BSET( IWORD, 73 
GOTO  81 

80  CALL  BCLR( IWORD, 73 

81  RETURN 
END 
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PROGRAM  READFILE 


C  PROGRAM  READFILE  READS  THE  DATA  FROM  THE  FILE 
C  CREATED  BY  PROGRAM  READTAPE.  THE  DATA  CAN  NOW  BE  PRO- 
C  CESSED  AS  THE  USER  WISHES.  THE  PICTURE  DATA  (PIXEL 
C  VALUES)  IS  STORED  IN  THE  ARRAY  MATRIX.  MATRIX(1,1)  IS 

C  THE  LEFT  REAR  OF  THE  IR  SCENE,  MATRIX ( 1 00 , 1 )  IS  THE 

C  RIGHT  REAR,  MATRIX  (100,1)  IS  THE  LEFT  FRONT,  AND, 

C  MATRIX( 1 00 , 1 00 )  IS  THE  RIGHT  FRONT 

C 

C 

COMMON  MATRIX( 10240) 

DIMENSION  BKGND(IO) ,AGENCY(10) , DATE( 3 ) ,PICCOM(35) 

+  DESC( 20) , INST (10) , SCENE  ( 35 ) ,C0ND(5) ,HIST(20) , 

INTEGER  TCODE , CMCODE , SITE , BEARING , BNCODE , BKGND, 

+  DATE, WEATHER, DESC, SCENE, COND, HIST, PICCOM, AGENCY 

INTEGER  DUMMY(64) ,DUM(7) ,0UTFILE(7) 

CALL  IOF( 1 . MAIN , OUTFILE , II  ,12, 13, MS, 14, 15, 16, 17) 

CALL  DFILWv "OUT" , I ) 

CALL  DFILW( "MATRIX", I) 

CALL  OPEN( 2 , "PICTURE" ,3,128,1) 

CALL  CFILW( "OUT" ,3,104,1) 

CALL  0PEN(3 , "0UT",3 ,1, 128) 

CALL  CFILW( "MATRIX", 3 ,40,1) 

CALL  OPEN( 4 , "MATRIX" , 3  » I) 

READ( 2 ,100)  TCODE, BKGND, CMCODE, BNCODE, AGENCY, SITE, 

+  DATE , WEATHER , DESC, INST, IVUANG , IDEPANG , IRANGE , SCENE , 

+  IRPM,COND,(HIST(I) ,1=1 ,10) ,(HIST(I) ,1=11 ,20) , 

+  BEARING, IWIND.IVTRAN1  , IVTRAN2 , IRTRAN1 , IRTRAN2 , 

+  IRC0DE1 , IRC0DE2 , IPYRA1 , IPYRA2 , IPYRHE1 ,IPYRHE2 

100  FORMAT ( 13 , 10A2 ,213 , 10A2 , 13 ,3A2 , A1 ,20 A2 , 10A2 ,16/216 , 

+  35A2,I6,5A2,10A2/10A2,2I6,10I3) 

READ( 2 , 1 20 )  TIME , CFSLOPE , CFINTCPT , PICCOM , HORRES , 

+  VERRES, HUMID, WIND, AIRTEMP , BARPRES , RAIN, SNOW, SOIL, 

+  VISABLE 

12U  FORMAT(/F10.3,2F20.15,35A2,F6.3/F6.3 ,2F1 0 .2 , FI  0 . 1 , 

+  FI  0 .0 , 4F1 0.2) 

READ( 2,160)  (MATRIX( J) , J=1 ,10000) 

160  FORMAT(//99(3(31I4/) ,714/) , 3( 31 14/ ) ,7 14 ) 

WRITE( 4,200)  MATRIX 
200  FORMAT(100(/"  ",3(2514/"  ")25I4)) 

CALL  WRBLK(4,0, MATRIX, 40, IE) 

IF(IE.NE.I)  TYPE  "READFILE  WRBLK4=",IE 

THE  PROGRAM  WILL  SKIP  OVER  THE  SECTION  THAT 
WRITES  DATA  TO  THE  FILE  CALLED  OUT. 

IIII=1 

IF(IIII. EQ.O)GOTO  666 
GOTO  333' 

C 

666  CONTINUE 
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WRITE( 3 ,110)  PICCOM, TCODE , SCENE , CMCODE 

110  F0RMAT("0",35A2,10X, "TARGET  CODE= " , 13/ "0" , 35 A2 , 1 0X, 
+  "COUNTERMEASURES  CODE=",I3) 

WRITE(3,111 )  BKGND,BNCODE,DESC, AGENCY 

111  F0RMAT("0", "BACKGROUND  -  ", 10A2 ,47X, "BANDCODEr  ”,I3 

+  /"0",20A2,40X, "AGENCY-”, 10A2) 

WRITEC 3 , 1 1 2 )  DATE , TIME , SITE , IVUANG , INST 

112  FORMAT ( ”0" , "DATE-" , 1 X , 3 A2 , 10X, "TIME-",  IX, FI 0 .3 , 

+  42X, "SITE-", 13, /"O", "VIEWING  ANGLE-", 14,"  DEG", 

+  5SX, "MEASURING  INST  ",10A2) 

WRITE( 3 ,113)  IDEPANG, IPYRA1 , IPYRA2 , BEARING, IPYRHE1  , 
+  IPYRHE2 

113'  FORMAT("0", "DEPRESSION  ANGLE-", 13,"  DEG",56X, 

+  "PYRANOMETER  READING-" , 1 X , 213 , "  W/M*»2"/"0", 

+  "BEARING-", 14 ,"  DEG" ,64X, "PYRHEL IOMETER  READING-", 

+  IX, 213,"  W/M**2" ) 

WRITE(3 ,114)  IRANGE, CFSLOPE , COND , CFINTCPT 

114  FORMATC'O",  "RANGE-", 13,"  METERS", 64X, 

+  "CONVERSION  FACTOR-SLOPE= " , F20 . 1 5/ 

+  "0", "PRIOR  VEHICLE  STATE-" ,5A2 ,50X, "CONVERSION 

+  FACTOR- INTERCEPT=" ,  F20.15) 

WRITEC 3,115)  HIST, SOIL, IWIND, WIND 

11b  FORMAT( "0" , 20A2 ,40X, "SOIL  MOISTURE  CONTENT? ",  FI  0 . 2/ 
+  "0" , "WIND  DIRECTION-", 14,"  DEG",7X, "WINDSPEED", 

+  F10.2,”  KNOTS") 

WRITE(3 ,116)  IVTRAN1 , IVTRAN2 , IRTRAN1 , IRTRAN2 , 

+  VISABLE 

11b  FORMATC'O",  "VISABLE  TRANSMISSION-"  ,213  ,  " 

+  "IR  TRANSMISSION-", 213, IX, "? "/ "0" , "VISABILITY" , 

+  FI  0 .2 , 1 X , "KM") 

WRITE(3 ,117)  HUMID, AIRTEMP,BARPRES, RAIN, SNOW 

117'.  FORMATCO",  "RELATIVE  HUMIDITY-" ,  FI  0 . 2 , 1  X,  "?"/ "0" , 

+  "AIR  TEMPERATURE-", IX, F10.1 ,"  DEG  CENTIGRADE "/ "0" , 
"BAROMETRIC  PRESSURE  ",  FI  0 . 0 , "  MILLIBARS"/ "0 » , 

+  "RAIN  RATE-" , FI  0 . 2  ,  "  MM  PER  H0UR"/"0", 

+  "SNOW  TEMPERATURE-", IX, FI 0.2, "DEG  CENTIGRADE") 
WRITE( 3 ,118)  HORRES , VERRES 

11b  FORMATC "0", "HORIZONTAL  RESOLUTION", 

+  1 X, F6 .3 , w  MRAD", 10X, "VERTICLE  RESOLUTION  ",F6.3, 

+  "  MRAD") 

WRITE( 3 ,122) 

122  FORMATC "0",45X, "A  VALUE  OF  999999  INDICATES  NO", 

+  "  DATA  AVAILABLE") 

DO  121  K=1 ,100 

WRITEC 3, 11 9)  K,CMATRIXC J,K) , J=1 ,100) 

121  CONTINUE 

119  FORMATC "0", "ROW", 14,/"  ",3(2515/"  "),25I5) 

333  CONTINUE 

WRITEC 10, 99) 0UTFILEC1) 

99  FORMATC"  ",S13, "created  by  READFILE") 

TYPE"MATRIX  CREATED  BY  READFILE" 

CALL  RESET 
END 
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PROGRAM  RMAT 


Program  RMAT  converts  a  100x100  TABILS  format 
picture  matrix  into  a  100x100  real  matrix 


DIMENSION  10(512) ,MATRIX( 1 0240) 

REAL  MA(256) 

COMMON  MA 

EQUIVALENCE(MA(1 ) ,10(1 )) 

INTEGER  OUTFILE(7) ,INFILE(7) 

CALL  IOF( 2, MAIN, INFILE , OUTFILE, II , 12 , MS , 13 , 14 , 15 , 16) 

2  FORMATC "  ",  SI  4  ,  "  created  by  RMAT") 

CALL  OPENC4, INFILE, 1  ,IE) 

IF(IE.NE.I)  TYPE  "RMAT  OPEN4=",IE 
CALL  DFILW( OUTFILE , IE) 

IF(IE.NE.I)  TYPE  "RMAT  DFILW=»,IE 
CALL  CFILW(0UTFILE,3,80,IE) 

IF(IE.NE.I)  TYPE  "RMAT  CFILW=«,IE 
CALL  OPEN(5, OUTFILE, 3, IE) 

IF(IE.NE.I)  TYPE  "RMAT  OPEN5=",IE 

910  FORMAT(100(/3(25I4/)25I4)) 

READ( 4,910)  MATRIX 

CALL  RDBLK(4,0, MATRIX, 40, IE) 

IF(IE.NE.I)  TYPE  "RMAT  RDBLK4=«,IE 

SUM=0 .0 

DO  11  J  =  1 ,10000 
SUM=SUM+MATRIX( J) 

11  CONTINUE 

SUM=SQRT( SUM) 

KK=0 

DO  100  1=1  ,40 
DO  200  J=1  ,256 
200  RA( J) =0 .0 

DO  400  J=1 ,256 
KK=KK+1 

400  MA( J) =MATRIX( KK)/SUM 

300  CALL  WRBLK(5,((I*2)-2) ,10, 2, IE) 

IF( IE. NE. 1 )TYPE  "RMAT  WRBLK=",IE 
100  CONTINUE 

WRITE( 10,2)  OUTFILE(I) 

CALL  RESET 
END 
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PROGRAM  PIX1 


Program  PIX1  converts  an  input  256x256  real 
number  file  into  a  VIDEO  file. 

Version  2 


INTEGER  0UTFILEC7) 

DIMENSION  10(2048) ,101(1024) ,102(256) ,INFILE(7) 
COMMON  RMAG( 1 024 ) 

EQUIVALENCE( RMAG( 1 ) , I0( 1 ) ) 

CALL  I0F(2, MAIN, INFILE, 0UTFILE, II  , 12 , MS , 13 , 14 , 15 , 16 ) 

9  FORMATC'  ",S13,"  created.") 

CALL  0PEN(3, INFILE, 1 , IE) 

IF(IE.NE.I)  TYPE" PI XI  0PEN3=",IE 
CALL  DFILW( OUTFILE , IE) 

IF(IE.NE.I)  TYPE"PIX1  DFILW=",IE 
CALL  CFILW( OUTFILE , 3 ,64 , IE) 

IF(IE.NE.1)TYPE" PIX1  CFILW=",IE 
CALL  0PEN( 4 , OUTFILE, 3 , IE) 

IF(IE.NE.I)  TYPE  "PIX1  0PEN4=",IE 

Find  min  and  max  values  of  infile, 

RMIN=1 .0E60 
RMAX=0 .0 
DO  2  1=0,63 

CALL  RDBLK( 3, (8*1), 10, 8, IE) 

IF(IE.NE.I)  TYPE "PIX1  RDBLK1 =" , IE 
DO  3  J=1 , 1 024 

IF(RMAG(  J) .GT.RMAX)  RMAX=RMAG( J) 

IF(RMAG( J) .LT.RMIN)  RMIN=RMAG( J) 

3  CONTINUE 

2  CONTINUE 

TYPE  "MIN= " , RMIN , "  MAX=",RMAX 

REWIND  3 

Convert  reals  to  gray  scale  integers  and  pack. 

ACCEPT  "PIX1  gray  scale  max?  ",RMAX1 
X  RMAX1 =350 

IF( RMAX. GT. RMAX1 )  RMAXsRMAXI 
ACCEPT  "PIX1  gray  scale  rain?  ",RMIN1 
X  RMIN1 =90 

IF( RMIN. LT.RM INI)  RMINsRMINI 
DO  5  1=0,63 

CALL  RDBLK(3,(I*8),IO,8,IE) 

IF(IE.NE.I)  TYPE  "PIX1  RDBLK2=" , IE 

6  4 


DO  4  J=1 ,1024 

IF(RMAG(J) .GT.RMAX)  RMAG( J) =RMAX 
IF(RMAG( J) .LT.RMIN)  RMAG( J) =  RMIN 
A=15 .0*( RMAG( J) -RMIN)/ ( RMAX-RMIN) 
I01( J)=IFIX(A) 

IF(I01(  1) .GT.15)  GOTO  9 00 

4  CONTINUE 
KK=0 

DO  6  K=1 ,256 
I02( K) =0 
DO  V.  J=1  ,4 
KK=KK+1 

I02(K)=ISHFT( 102 (K) ,4) 

7  102 (K) =102 ( K) +101 ( KK) 

6  CONTINUE 

CALL  WRBLfC(  4 , 1 , 102 , 1  ,  IE) 
IF(IE.NE.I)  TYPE "PIX1  WRBLK=",IE 

5  CONTINUE 

WRITE( 10,9)  OUTFILE(I) 

CALL  RESET 
GOTO  901 

900  TYPE  "PIX1  gray  scale  error” 

901  CONTINUE 

END 


X 
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PROGRAM  CMAT 


Program  CMAT  converts  a  100x100  TABILS  format 
picture  matrix  into  a  256x256  complex 
array 


DIMENSION  10(1024), MATRIX( 10240) 

COMPLEX  MA( 256 ) 

COMMON  MA 

EQUIVALENCE(MA(1) ,10(1)) 

INTEGER  0UTFILE(7) ,INFILE(7) 

CALL  I0F( 2, MAIN, INFILE , OUTFILE , II , 12 , MS , 13 , 14 , 15 , 16 ) 

2  FORMATC  ",S14,"  created  by  CMAT") 

CALL  0PEN(4, INFILE, 1  , IE) 

IF(IE.NE.I)  TYPE  "CMAT  OPEN4=",IE 
CALL  DFILVK  OUTFILE ,  IE) 

IF(IE.NE.I)  TYPE  "CMAT  DFILW=",IE 
CALL  CFILW(OUTFILE,3,1024,IE) 

IF(IE.NE.I)  TYPE  "CMAT  CFILW=",IE 
CALL  0PEN(5, OUTFILE, 3, IE) 

IF(IE.NE.I)  TYPE  "CMAT  0PEN5=",IE 

910  F0RMAT(100(/3(25I4/)25I4)) 

READ(4 ,910)  MATRIX 

CALL  RDBLK( 4,0, MATRIX, 40, IE) 

IF(IE.NE.I)  TYPE  "CMAT  RDBLK4=" , IE 

KK=0 

DO  100  1=1,256 
DO  200  J=1 ,25b 
20U  MAC J)=(0. 0,0.0) 

IF(I.LE.78.0R.I.GT.178)  GOTO  300 

DO  400  K=79  » 17  8 

KK=KK+l; 

IFCMATRIX(KK) .GT.1023)  MATRIX(KK)=0 
A=MATR1X( KK) 

400  MA( K) =CMPLX( A, 0 .0) 

300  CALL  WRBLKC 5,(CI*4)-4) ,10, 4, IE) 

IF(IE.NE.1)TYPE  "CMAT  WRBLK=",IE 
100  CONTINUE 

WRITEC 10,2)  OUTFILE(I) 

CALL  RESET 
END 
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PROGRAM  CREAD 


PROGRAM  CREAD  READS  VALUES  FROM  CMATRIX,  AND 
CONVERTS  THEM  INTO  A  MAGNITUDE  AND  PHASE  FORMAT 
MAGNITUDES  ARE  WRITTEN  TO  A  FILE  RMAG. 

PHASE  ANGLES  ARE  WRITTEN  TO  A  FILE  CALLED  ANG. 


DIMENSION  10(1024) ,RMAG( 256) ,ANG(256) ,101(512)  , 
102(512) 

COMPLEX  MA(256) 

COMMON  MA, RMAG, ANG 

EQUIVALENCE (MA( 1 ) ,I0(  1 ) ) 

EQUIVALENCE( 101(1), RMAG(  1  ) ) 

EQUIVALENCE( 102(1) ,ANG(1)) 

DIMENSION  INFILE( 7 ) 

CALL  IOF( 1 ,MAIN, INFILE, II , 12 , 13 , MS , 14 , 15 , 16 , 17) 


CALL  DFILW("RMAG",I) 

IF(I.NE.I)  TYPE  "CREAD  DFILW  RMAG=",I 
CALL  DFILW( "ANG" , I) 

IF(I.NE.I)  TYPE"CREAD  DFILW  ANG=",I 
CALL  CFILW( "RMAG" ,3,512,1) 

IF(I.NE.I)  TYPE  "CREAD  CFILW  RMAG=",I 
CALL  CFILW( "ANG", 3 ,512,1) 

IF(I.NE.I)  TYPE"CREAD  CFILW  ANG=",I 
CALL  0PEN(3, INFILE, 1 , I) 

IF(I.NE.I) TYPE "CREAD  0PEN3=",I 
CALL  OPEN( 4 , "RMAG" , 3,1) 

IF(I.NE.I)  TYPE"CREAD  0PEN4=",I 
CALL  OPEN(5,"ANG",3,I) 

IF(I.NE.I)  TYPE"CREAD  0PEN5=",I 

DO  1  1=0,255 

CALL  RDBLK(3,(I*4) ,10, 4, IE) 
IF(IE.NE.I)  TYPE*  CREAD  RDBLK=",IE,I 
DO  2  J=1  ,256 
RMAG( J ) =CABS( MA( J) ) 

IF(RMAG( J) .EQ.O.O)  GOTO  44 
X=REAL( MA( J) ) 

Y=AIMAG(MA( J) ) 

ANG( J ) =ATAN2( Y, X) 

GOTO  2 
ANG( J) =0 .0 
CONTINUE 

CALL  WRBLK( 4, (I *2) ,101 ,2, IE) 
IF(IE.NE.I)  TYPE"CREAD  WRBLK4=",IE 
CALL  WRBLK(5,(I»2) ,102, 2, IE) 
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IF(IE.NE.I)  TYPE  "CREAD  WRBLK5=",IE 
1  CONTINUE 
CALL  RESET 

STOP  Magnitude  and  phase  files  created.  RMAG,  ANG. 
END 


t 
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PROGRAM  HOREV 


Program  HOREV  takes  the  magnitude  file  from 
the  template  image  and  combines  it  with  the  phase 
file  from  the  image  and  combines  them  to  make  a 
modified  256x256  complex  image  file. 


COMPLEX  MA( 256 ) 

DIMENSION  RMAG(256; ,ANG(256) 

COMMON  101(1024) ,102(512) ,103(512) 

INTEGER  INFILE1  (7) ,INFILE2(7) ,0UTFILE(7) 
EQUIVALENCE(IOKI)  ,MA(1)) 

EQUIVALENCE( I02( 1 ) , RMAG( 1 ) ) 

EQUIVALENCE( 103(1) , ANG( 1 ) ) 

C 

CALL  I0F(3,MAIN,INFILE1  , INFILE2 , OUTFILE , I , M , I , I , I , I) 
C 

2  FORMATC  ",S13, "Created  by  HOREV") 

C 

CALL  DFILW( OUTFILE, IE) 

IF(IE.NE.I)  TYPE"HOREV  DFILW=",IE 
CALL  CFILfcK OUTFILE, 3 , 1024 , IE) 

IF(IE.NE.I)  TYPE"HOREV  CFILW=",IE 
CALL  0PEN(3, OUTFILE, 3 ,IE) 

IF(IE.NE.I)  TYPE"HOREV  0PEN3=",IE 
CALL  0PEN(4',INFILE1  ,1  ,IE) 

IF(IE.NE.I)  TYPE"HOREV  0PEN4=",IE 
CALL  OPEN( 5 , INFILE2 , 1 ,IE) 

IF(IE.NE.I)  TYPE"HOREV  0PEN5=",IE 
C 

DO  4  1=0,255 

CALL  RDBLK(4,(I«2) ,102, 2, IE) 

IF(IE.NE.I)  TYPE "HOREV  RDBLK4= " , IE , I 
CALL  RDBLK( 5 , ( 1*2) , 103 ,2 , IE) 

IF(IE.NE.I)  TYPE"HOREV  RDBLK5=", IE, I 

DO  3  J  =  1  ,256 

X=RMAG( J) *COS( ANG( J) ) 

Y=RMAG( J) *SIN( ANG( J) ) 

MA( J) =CMPLX( X, Y) 

3  CONTINUE 

CALL  WRBLK(3,(I*4) ,101 ,4, IE) 

IF(IE.NE.I)  TYPE"HOREV  WRBLK=",IE,I 

4  CONTINUE 
CALL  RESET 

WRITE( 10,2)  OUTFILE(I) 

END 
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PROGRAM  FILTER 


Program  filter  calls  a  complex  256x256  file 
and  filters  it  by  a  function  F( ROW, COL)  specified 
by  the  user.  ROW  equals  the  absolute  value  of 
the  distance( number  of  rows)  from  the  129th  row  of 
the  array.  COL  equals  the  absolute  value  of  the 
distance( number  of  columns)  from  the  129th  column 
Thus  R0W=0,  COL=0  specifies  the  element  A(129,129) 
of  array  A.  Also  R0W=1 ,  COL=0  specifies  the  ele¬ 
ments  A(128,129)  and  A(130,129)  of  the  array  A. 

To  change  F( ROW, COL)  the  user  must  change  the 
function  of  ROW  and  COL  in  this  program  and 
recompile  and  reload  the  program. 


COMPLEX  MA( 256 ) 

COMMON  10(1024) 

INTEGER  INFILEC7) ,0UTFILE(7) 

EQUIVALENCE( I0( 1 ) ,MA( 1 ) ) 

CALL  I0F( 2 . MAIN , INFILE , OUTFILE , II , 12 , MS , 13 , 14 , 15 , 16 ) 
CALL  DFILW( OUTFILE, IE) 

IF(IE.NE.1 .AND. IE. NE. 13) GOTO  12 
CALL  CFILW( OUTFILE, 3, 1024, IE) 

IF(IE.NE.I)  GOTO  13 
CALL  0PEN(0, OUTFILE, 3, IE) 

IF( IE.NE. 1 )G0T0  14 
CALL  0PEN(1 , INFILE, 1 , IE) 

IF(IE.NE.I)  GOTO  15 
DO  1  K=0 ,255 

CALL  RDBLKC 1 ,K*4 ,10,4 , IE) 

IF( IE. NE. 1 ) GOTO  10 
DO  2  J=0 ,255 
IF( J . GT. 1 28) GOTO  3 
ROW= 1 28- J 
GOTO  4 

3  ROW= J-128 

4  CONTINUE 
IFCK.GT.128) GOTO  5 
COL= 1 28-K 

GOTO  6 

5  C0L=K-128 

6  CONTINUE 

xxx**xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

XXXXAttXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

THE  FILTER  FUNCTION  SHOULD  BE  CONTAINED  IN  THIS 

SECTION 

X**«*ftXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXXX 

MA(J+1)=MA(J+1)*SQRT( ROW**2+COL**2)  ;  F( ROW, COL) 


cxx»***xxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxxx 

(;XXX«*t**tXX«XXXXXXXtXiXXXXXXXXXXXXXXXXXXXXXX*XXXXXXXXXXXXXX 

2t  CONTINUE 
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CALL  WRBLK( 0 , K*4 , 10 , 4 , IE) 

IF( IE. ME. 1 )G0T0  11 
1  CONTINUE 

WRITEOO ,9)  OUTFILE(I) 

9  FORMAT(  "  ,  S 1 3  ,  "created  by  FILTER”) 

GOTO  999 


10 

STOP”FILTER 

RDBLK=”,IE,K 

11 

STOP”FILTER 

WRBLK= ” , IE , K 

12 

STOP”FILTER 

DFILW=  n , IE 

13 

STOP”FILTER 

CFILW=M, IE 

14 

STOPnFILTER 

OPENOs", IE 

15 

STOP”FILTER 

OP EN1s",IE 

999. 

CONTINUE 

END 
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PROGRAM  CORRELATE 


Program  CORRELATE  multiplys  two  256x256  com¬ 
plex  files  according  to  the  formula; 

INF1LE1 (  J, K) *CONJ(  INFILE2C  J , K) ) =OUTFILE( J , K) . 
The  resulting  array  is  stored  in  file  OUTFILE. 


DIMENSION  INFILEK7)  ,INFILE2(7) 

INTEGER  0UTFILE(7) 

COMPLEX  MAI  (256) ,MA2(256) ,MA3( 256) 

COMMON  101(1024), 102(1024) ,103(1024) 

EQUIVALENCE ( MAI ( 1 ) ,101  (  1 ) ) 

EQUIVALENCE( MA2( 1 ) ,102(1)) 

EQUIVALENCE(MA3(1) ,103(1)) 

C 

CALL  I0F(3 ,MAIN, INFILE 1 , INFILE2 , OUTFILE , II ,MS,I2, 
+13*14,15) 

C 

2  FORMATC  ,f ,  S 1  3 "Created  by  CORRELATE”) 

C 

CALL  DFILW( OUTFILE, IE) 

IF(IE.NE.I)  TYPE "CORRELATE  DFILW=",IE 
CALL  CFILW( OUTFILE , 3 , 1 024 , IE) 

IF( IE. NE . 1 ) TYPE" CORRELATE  CFILW=" , IE 
CALL  OPEN(1 , OUTFILE, 2, IE) 

IF( IE. NE.1) TYPE "CORRELATE  0PEN1=",IE 
CALL  OPEN( 2 , INFILE1 ,2, IE) 

IF( IE. NE. 1 ) TYPE "CORRELATE  0PEN2=", IE 
CALL  OPEN ( 3 , INFILE2 ,2 , IE) 

IF( IE. NE. 1) TYPE "CORRELATE  0PEN3=",IE 
C 

DO  4  1=0,255 

CALL  RDBLK( 2, (I *4) ,102, 4, IE) 

IF(IE.NE.I)  TYPE"CORRELATE  RDBLK2= " . IE , I 
CALL  RDBLK(3,(I*4),I03,4,IE) 

IF( IE. NE. 1 ) TYPE "CORRELATE  RDBLK3 = " , IE , I 
DO  3  J=  1  ,256 

MAI ( J) =MA2( J ) *CONJG( MA3( J) ) 

3  CONTINUE 

CALL  WRBLK( 1,(1 *4), 10 1,4, IE) 

IF( IE . NE . 1 ) TYPE "CORRELATE  WRBLK1 =" , IE , I 

4  CONTINUE 

WRITE( 10,2)  OUTFILE(I) 

END 
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PROGRAM  RBIG 


Program  RMAT  takes  9  scenes  created  by  RMAT  and 
combines  them  into  one  large  256x256  real  scene 
array. 


DIMENSION  IFILEC7) 

CALL  IOF( 1 , MAIN, IFILE,  12, 13, 14, MS,  15 ,16,17,18) 
CALL  DFILW(IFILE,IE) 

IFQE.NE.1  .AND.  IE.  NE.  13)  TYPE  "BIG  DFILW=",IE 
CALL  CFILW( IFILE, 3, 512, IE) 

IF( IE . NE . 1 ) TYPE "BIG  CFILW=",IE 
CALL  OPEN ( 1 , "MATR1X1 ",  1 , IE) 

IF( I£ . NE. 1 ) TYPE"BIG  OPEN1=",IE 
CALL  OPEN( 2 , "MATR1X2" , 1 ,IE) 

IF(  IE.  NE.  1  )TYPE''BIG  0PEN2=",IE 
CALL  OPEN(3 , "MATR1X3" , 1 ,IE) 

IF( IE. NE. 1 )TYPE"BIG  OPEN3=",IE 
CALL  OPEN (4, IFILE, 3, IE) 

IF( IE. NE. 1 ) TYPE"BIG  OPEN4=",IE 
CALL  READ3U  ,2, 3, 4,1) 

CALL  RESET 

CALL  OPEN ( 1 , "MATR1X4" ,1 , IE) 

IF( IE. NE. 1 ) TYPE"BIG  0PEN5=",IE 
CALL  OPEN( 2 , "MATR1X5 " , 1 ,IE) 

IF(  IE.  NE.  1 )  TYPE"BIG  0PEN6=",IE 
CALL  OPEN(3,"MATRIX6",1 ,IE) 

IFCIE.NE.1 )TYPE"BIG  OPEN7=",IE 
CALL  OPEN( 4, IFILE, 3, IE) 
IF(IE.NE.1)TYPE"OPENJJ2=",IE 
CALL  READ3(1 ,2, 3,4, 2) 

CALL  RESET 

CALL  OPEN(1 , "MATRIX7 " , 1 ,IE) 

IF( IE. NE. 1 ) TYPE"BIG  OPEN8=",IE 
CALL  OPEN( 2 , "MATR1X8" , 1 ,IE) 
IF(IE.NE.1)TYPE"OPEN9=",IE 
CALL  OPEN(3,"MATRIX9",1 ,IE) 

IFCIE.NE.1 ) TYPE "BIG  OPEN13=",IE 
CALL  OPEN(4,IFILE,3,IE) 

IF(IE.NE.I) TYPE "OPEN  JJ=3",IE 
CALL  READ3C1 ,2, 3, 4, 3) 

WRITE( 10,44)  IFILE(I) 

44  FORMAT( "  ",  SI  3 , "created  by  RBIG") 

END 
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PROGRAM  READ3 


C  Subroutine  READ3  reads  in  3  picture  arrays  created 
C  by  RMAT  and  writes  them  to  the  9x9  picture  array  being 
C  created  by  RBIG.  That  is  the  9x9  large  scene  is  created 
C  one  "row"  of  pictures  at  a  time. 

DIMENSION  RA(256) 

INTEGER  MA4(512) 

COMMON  MA4 

EQUIVALENCE ( RA( 1 ) ,MA4( 1 ) ) 

REAL  MA1( 100) ,MA3( 100) ,MA2(100) 

SUBROUTINE  READ3( CHI , CH2 , CH3 , CH4 , J J ) 

INTEGER  CHI ,CH2,CH3,CH4 
DO  4  J=1 ,80 
READ  BINARY(CHI)  MAI 

READ  BINARY( CH2)  MA2 

READ  BINARY(CH3)  MA3 

IF( J.LE.5)  GOTO  4 
DO  1  1=1  ,78 

1  RA( I) =MA1 (I) 

DO  2  1=79,172 

2  RA( I) =MA2( 1-78) 

DO  3  1=173,256 

3  RA( I) =MA3( 1-172) 

L=2*( J-5*( J J-1 ) -1 ) ;  increment 

KK=L+160*( JJ-1 ) 

CALL  WRBLK(CH4,KK,MA4,2,IE) 

IF(IE.NE.1)TYPE«3READ  WRBLK=",IE,J, JJ 

4  CONTINUE 
RETURN 
END 
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PROGRAM  IOF 


SUBROUTINE  IOF< N , MAIN , FI , F2 , F3 , F4 , MS , SI ,S2,S3,S4) 

Written  by  Lt.  Simmons  31  Aug  1981 

This  FORTRAN  5  subroutine  will  read  from  the  file 
COM. CM  (FCOM.CM  in  the  foreground)  the  program  name, 
any  global  switches,  and  up  to  four  local  file  names 
and  corresponding  switches. 

Calling  arguments: 

N  is  the  number  of  local  files  and  switches  to  be 
read  from  (F)COM.CM.  N  must  be  1,  2,  3,  or  4. 

MAIN  is  an  array  for  the  main  program  file  name. 

FI,  F2,  F3 ,  and  F4  are  the  four  variables  to  return 
the  local  file  names. 

MS  is  a  two-word  integer  array  that  holds  any  global 
switches . 

SI,  S2,  S3,  and  S4  are  two-word  integer  arrays  that 
hold  the  local  switches  corresponding  to  FI  through 
F4  respectively. 

Dimension  the  arrays. 

DIMENSION  MAIN(7) ,MS(2) 

INTEGER  FI (7) ,F2(7) , F3(7) ,F4(7) ,S1 (2) 

+  ,S2(2) , S3 (  2 ) ,S4(2) 

Check  the  bounds  on  N. 

IF( N. LT. 1 . OR . N . GT. 4 ) STOP  "N  out  of  bounds  in  IOF" 

Process  the  data  in  COM. CM  (or  FCOM.CM). 

CALL  GROUND(I)  ;Find  out  which  ground  I  am  in 
IF( I. EQ. 0) OPEN  0, "COM. CM"  ;0pen  ch.  0  to  COM. CM 

IF(I.EQ. 1 )0PEN  0, "FCOM.CM"  ;0pen  ch.  0  to  FCOM.CM 

CALL  C0MARG( 0 , MAIN , MS , IER )  ;Read  from  COM. CM 

IF( IER. NE.  1 ) TYPE"  COMARG  error:", IER 
WRITE( 10 , 1 )MAIN( 1 )  ;Type  prog,  name 

1  F0RMAT(  ’  Program  * , SI  3 running. ' ) 

CALL  C0MARG(0,F1 ,S1 ,JER)  ;Read  from  COM. CM 

IF( JER. NE.  1 ) TYPE"  COMARG  error  (F1):",JER 
IF( N. EQ. 1 ) GO  TO  2  ;Tesl  N 

CALL  C0MARG( 0 ,F2 , S2 ,KER)  ;Read  from  COM. CM 

IF(KER. NE. 1 ) TYPE"  COMARG  error  (F2):",KER 
IF( N.EQ.2)G0  TO  2  ;Test  N 
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CALL  COMARG(0,F3,S3,LER)  ;Read  from  COM. CM 

IF( LER. NE. 1 )TYPE"  COMARG  error  (F3):",LER 
IF( N. EQ. 3 ) GO  TO  2  ;Test  N 

CALL  C0MARG(0,F4,S4 ,MER)  ;Read  from  COM. CM 

IF(MER.NE. 1 )TYPEM  COMARG  error  (F4):",MER 
2  CLOSE  0 
RETURN 
END 
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PROGRAM  UNPACK 


PROGRAM  UNPACK:  Program  unpacks  a  64X256  Video 
picture  into  a  256Pixel  by  256' array  located 
in  file  UFILE. 

SUBROUTINE  UNPACK( NAME , IFILE) 

DIMENSION  NAME(7) 

INTEGER  IPICT( 256 ) ,NPICT(1024) ,A,B 
CALL  CFILW  (IFILE, 2, IER) 

CALL  OPEN  ( 1 , NAME,0 , IER) 

CALL  OPEN  (2, IFILE, 0, IER) 

DO  1  1=0,63 
L=I*4 

CALL  RDBLK  ( 1 , I , IPICT , 1 , IER) 

DO  2  J  =  1  ,256 
M=  J  *4 

N=IPICT( J) 

DO  3  K=1  ,4 
A=15 .AND. N 
NPICT(M)=A 
N=ISHFT( N,-4) 

M=M-1 

3  CONTINUE 

2  CONTINUE 

CALL  WRBLK  ( 2 ,L, NPICT,4 , IER) 

1  CONTINUE 

CALL  RESET 
RETURN 
END 
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PROGRAM  REPACK 


V 


PROGRAM  REPACK:  Program  takes  a  normalized 
integer  file  (NIFILE)  and  packs  it  into 
a  file  (PICT)  usable  by  Video. 

SUBROUTINE  REPACK( IFILE) 

DIMENSION  IPICTC1024) ,NPICT(256) 

INTEGER  A,B,C 

CALL  CFILW  ( "PICT" ,2 ,IER) 

CALL  OPEN  (1 ,IFILEfO,IER) 

CALL  OPEN  (2,"PICT",0,IER) 

DO  1  1=0,63 
M=I«4 

CALL  RDBLKC1 , M, IPICT , 4 , IER ) 

DO  2  J=1 ,1024,4 
L=J 

N=( J+33/4 
A=0 

DO  3  K=1  ,4 

A=ISHFT( A , 4 ) 

A=IPICT(L)+A 
L=L+1! 

,  3  CONTINUE 

NPICT(N)=A 

2  CONTINUE 

CALL  WRBLK( 2 , I , NPICT , 1 ,IER) 

1  CONTINUE 

CALL  RESET 
RETURN 
END 
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PROGRAM  REDNOSE 


PROGRAM  REDNOSE:  This  program  performs 

a  user  specified  averaging  on  previously 
created  video  files  of  identical 
images.  Inputs  are  specified  for  each 
iteration  by  the  user.  The  final 
image  is  output  to  file  VAVG. 


DIMENSION  NAMEC7) 

INTEGER  FPICT(256) , SPICT( 256 ) , NPICT( 256 ) 

C 

K=0 

ACCEPT"  Name  of  first  file?  » 

READ( 11 ,15)  NAME(1) 

15  FORMAT  (SI  3) 

CALL  DUNPACK  ( NAME , "VPICT1 " ) 

TYPE"  File  is  unpacked  and  in  VPICT1," 

C 

ACCEPT"  Name  of  second  file?  " 

READC 11,15)  NAME(1) 

CALL  DUNPACK  ( NAME , "VPICT2") 

TYPE"  File  is  upacked  and  in  VPICT2.  " 

CALL  CFILW  ( "VAVG" , 2 ,  IER) 

IF(IER.NE. 1 )TYPE"File  create  error.", IER 
CALL  OPEN  (1 , "VAVG", 2, IER) 

IFCIER.NE.1 ) TYPE "Channel  1  open  error  1.",IER 
C 

CALL  OPEN  ( 2 , "VPICT1 " , 2 , IER ) 

IF( IER. NE. 1 ) TYPE "Channel  2  open  error  2.", IER 
CALL  OPEN  ( 3 , "VPICT2" ,2 , IER } 

IF(IER.NE.I) TYPE "Channel  3  open  error  3.", IER 
DO  1  1=0,255 

CALL  RDBLK  (2,I,FPICT, 1 ,IER) 

IF( IER. NE.1) TYPE "First  RDBLK  error.", IER 
CALL  RDBLK  (3 , I , SPICT, 1 , IER ) 

IF(IER.NE. 1 ) TYPE "Second  RDBLK  error.", IER 
DO  2  J= 1,256 

NPICT( J)=(FPICT( J)+SPICT( J) )/2 

2  CONTINUE 

CALL  WRBLK  ( 1  , 1 , NPICT , 1 , IER) 

IF(IER. NE. 1 ) TYPE "First  WRBLK  error.", IER 
1  CONTINUE 

K=K+2 

CALL  RESET 
C 

3  ACCEPT"  Do  you  wish  to  continue 
READ  (11,16)  IYN1 

1b  FORMAT  (SI) 

IF  ( IYN1 .EQ. 19968)  GO  TO  9 
IF  (IYN1 .NE. 22784)  GO  TO  7 
GO  TO  8 


(YES/ NO)?  " 


;N0  TO  STOP 
,'NOT  YES 
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7 

C 

C 

8 


17 


1  u 


c 

c 

9 

12 


C 

14 

13 


TYPE"  Input  error!!!  Try  again. 
GO  TO  3 


ACCEPT"  Name  of  next  file?  " 

READ( 11 ,15)  NAME(1) 

CALL  DUNPACK( NAME , "VPICTN" ) 

TYPE"  File  is  unpacked  in  VPICTN." 

CALL  OPEN  (1 ,"VAVG",2,IER) 

IF( IER. NE. 1 ) TYPE "Second  OPEN  channel  1  error. ",IER 
CALL  OPEN  (2, "VPICTN", 2, IER) 

IF( IER. NE. 1 ) TYPE"2nd  Channel  2  error.", IER 
CALL  CFILW  ( "VAVGN" ,2 , IER) 

IF( IER. NE. 1 ) TYPE  "  2nd  create  error.", IER 
CALL  OPEN  (3, "VAVGN", 2, IER) 

IF( IER. NE. 1 )TYPE"  2nd  open  channel  3  error.", IER 
DO  10  1=0,255 

CALL  RDBLK  ( 1 , I , FPICT , 1 , IER ) 

IF(IER.NE.I) TYPE "Third  RDBLK  error.", IER 
CALL  RDBLK  ( 2 , I , SPICT , 1 , IER) 

IFQER.NE.1  )TYPE"4th  RDBLK  error.", IER 
DO  17  J=1,256 

NPICT(J)=(FPICT(J) +SPICT( J ) )/2 
CONTINUE 

CALL  WRBLK  ( 3 , I , NPICT , 1 , IER) 

IF( IER. NE. 1 ) TYPE"2nd  WRBLK  error.", IER 
CONTINUE 
K=K+1 ; 

CALL  RESET 

CALL  DFILW( "VAVG",IER) 

IF( IER. NE. 1 )TYPE"VAVG  Delete  error.", IER 
RENAME  "VAVGN", "VAVG" 

IF( IER. NE. 1 ) TYPE "Rename  failed.", IER 
GO  TO  3* 


CALL  DREPACK( "VAVG") 

ACCEPT"  Do  you  wish  to  rename  the  averaged  file(YES/NO)?  ” 
READ  (11,16)  IYN2 

IFQYN2.EQ. 22784)  GO  TO  13'  ;YES 

IF(IYN2.NE. 19968)  GO  TO  14  ;NO 

TYPE"  You  have  averaged  ",K,"  pictures.  " 

TYPE"  Averaged  file  is  in  file  PICT.  " 

GO  TO  18 

TYPE"Input  error.  Try  again.  " 

GO  TO  12 

ACCEPT"Desired  outputfile  name?  " 

READ  (11,15)  NAME(1) 

RENAME  "PICT", NAME 
CALL  DFILW( "VAVG", IER) 

CALL  DFILW( "VPICT1 " ,IER) 

CALL  DFILW( "VPICT2" , IER) 

CALL  DFILW( "VPICTN", IER) 

STOP 
END 
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PROGRAM  CLEAN 


PROGRAM  CLEAN;  REVISION  1:  Program 
finishes  job  started  by  THRES. 

DIMENSION  10(256) ,NAME(7) 

INTEGER  A,B 

ACCEPT"  Name  of  input  file?  " 

READ( 11 , 30 )  NAME ( 1 ) 

30  FORMATCS13) 

CALL  DUNPACK( NAME , "TPICT") 

CALL  OPEN  (1  "TPICT", 1 ,IER) 

IF(IER.NE.I) TYPE "OPEN 1  ERR0R=  ",IER 
CALL  CFILW  ( ”PUP",2,IER) 

IFCIER.NE. 1 ) TYPE "CREATE  ERR0R=  ",IER 
CALL  OPEN  (2, "PUP", 3 ,IER) 

IF( IER. NE . 1 ) TYPE"0PEN2  ERR0R=  ",IER 
ACCEPT"  Threshold  value?  ",NUMi: 

DO  1  1=0,255 

CALL  RDBLK  ( 1 , I , 10 , 1 , IER ) 

IF(IER.NE.1 )TYPE"RDBLK  ERRORr  ",IER 
DO  2  J=1 ,25b 
A=IO( J) 

IF(  A.  LT.  NUM1 )  GO  TO  3 
I0( J) =15 

2  CONTINUE 
GO  TO  5 

3  DO  4  K= 1  ,25b 
L=257-K 
B=IO(L) 

IF( B. LT. NUM1 ) GO  TO  5 
I0( L) =1 5 

H  CONTINUE 

5  CALL  WRBLK  ( 2 , I , 10 , 1 , IER) 

IF( IER . NE. 1 )TYPE"WRBLK  ERRORr  ",IER 
1  CONTINUE 

CALL  RESET 
CALL  TREPACK( "PUP") 

ACCEPT"  Name  of  output  file.  " 

READ( 11 ,30)  NAME(1) 

RENAME  "TPICTN" , NAME 
CALL  DFILW  ( "TPICT" , IER) 

CALL  DFILW  ("PUP", IER) 

STOP 

END 


( 
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PROGRAM  COMBINE 


PROGRAM  COMBINE,  VERSION  2:  Takes  a  scene 
consisting  of  any  single  grey  level  and 
combines  it  with  a  template. 

DIMENSION  NAME(7) ,10(256) ,NO(256) ,LO(256) 
INTEGER  A, B 

ACCEPT"  Name  of  template  file?  " 

READ( 11,30)  NAME(1) 

30  FORMAT( SI  3) 

CALL  DUNPACK( NAME, "TEMP") 

ACCEPT"  Name, of  reversed  PIMT  file?  ” 
READ( 1 1 ,30)  NAME(1) 

CALL  DUNPACKC NAME, "IMAGE") 

CALL  0PEN(1 ,"TEMP",1 ,IER) 

CALL  OPEN  ( 2 , "IMAGE" ,1 , IER) 

CALL  CFILW  ( "NEW", 2, IER) 

CALL  OPEN  (3 , "NEW", 3 ,IER) 

DO  1  1=0,255 

CALL  RDBLK  ( 1 , I , 10 , 1 , IER) 

CALL  RDBLK  ( 2 , I , NO , 1 , IER) 

DO  2  J=1 ,256 
A=IO( J) 

IF( A. EQ. 1 5 )G0  TO  3 
L0( J ) =N0( J) 

GO  TO  2 

3  L0( J) =A 

2  CONTINUE 

CALL  WRBLK  ( 3 , I , LO , 1 , IER) 

1  CONTINUE 

CALL  RESET 
CALL  TREPACK( "NEW") 

ACCEPT"  Name  of  output  file?  " 

READ( 11 ,30)  NAME(1) 

RENAME  "TPICTN" , NAME 
CALL  DFILW  ("TEMP", IER) 

CALL  DFILW  ( "IMAGE", IER) 

CALL  DFILW  ("NEW", IER) 

STOP 

END 
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PROGRAM  MOVE 


PROGRAM  MOVE:  This  program  moves  a 

template  to  a  new  location  by  way  of  a  2-D 
shift.  It  is  assumed  that  the  location  of 
the  template  is  known. 

DIMENSION  NAME(7) ,10(256) ,NO(256) 

INTEGER  A,B,C,D 

ACCEPT"  Name  of  template  to  be  moved?  " 
READ( 11 ,30)  NAME(1) 

30  FORMATC  SI  3 ) 

CALL  DUNPACK( NAME, "TPICT") 

CALL  OPEN  ( 1 , "TPICT", 1 ,IER) 

IF(IER.NE.1)TYPE"0PEN1  ERRORr  ",IER 
CALL  CFILW  ( "TEM",2,IER) 

IF(IER.NE. 1 )TYPE"CREATE  ERR0R=  ",IER 
CALL  OPEN  ( 2 , "TEM " , 3 , IER ) 

IFCIER.NE.1 )TYPE"0PEN2  ERR0R=  ",IER 
ACCEPT"  Enter  old  top  row  number.  ",NUM1 
ACCEPT"  Old  bottem  row  number?  ",NUM2 
ACCEPT"  Old  right  column  number?  ",NUM3 
ACCEPT"  Old  left  column  number?  ",NUM4 
ACCEPT"  New  top  row?  ",NUM5 
ACCEPT"  New  left  column  number?  ",NUM6 
i  ACCEPT"  Background  grey  level?  ",NUM7 

A=NUM1 -NUM5 
D=NUM3-NUM4 
B=NUM4-NUM6 
C=NUM2-NUM1 
NBR=NUM5+C 
NRC=NUM6+D 
DO  1  1=0 , NUM5 
DO  2  J=1 ,256 
NO( J)=NUM7 
2  CONTINUE 

CALL  WRBLK  ( 2 , 1 , NO , 1 , IER) 

IF(IER.NE.I) TYPE"WRBLK1  ERR0R=  ",IER 
1 ;  CONTINUE 

DO  3  I=NUM1 , NUM2 

CALL  RDBLK  ( 1 , I , 10 , 1 , IER) 

IF( IER. NE. 1 )TYPE"RDBLK  ERR0R=  ",IER 
M=I-A 

DO  4  J=1 ,NUM6 
N0( J)=NUM7 

4  CONTINUE 

DO  5  J=NRC,256 
NO( J)=NUM7 

5  CONTINUE 

DO  6  J=NUM4 , NUM3 
f  L= J-B 

N0( L) =IO( J) 

6  CONTINUE 
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CALL  WRBLK  ( 2  ,M, NO, 1 , IER) 

IFCIER.NE.1 )TYPE"WRBLK2  ERROR=  ",IER 
CONTINUE 
DO  7  I=NBR,255 
DO  8  J=1 ,256 
NO( J)=NUM7 
CONTINUE 

CALL  WRBLK  ( 2 , I , NO , 1 , IER) 
IF(IER.NE.1)TYPE"WRBLK3  ERROR=  »,IER 
CONTINUE 
CALL  RESET 
CALL  TREPACK( "TEM") 

CALL  DFILW( "TEM",IER) 

CALL  DFILW( "TPICT" , IER ) 

ACCEPT"  Desired  name  for  output  file?  " 
READ( 11 ,30)  NAME(1) 

RENAME  "TPICTN" , NAME 

STOP 

END 
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PROGRAM  EXAMINE 


PROGRAM  EXAM:  This  program  asks  for  a 

threshold  value.  If  the  pixel  is  less  than 
the  specified  threshold,  it  is  put  in 
the  new  image.  If  not,  the  pixel  is  set 
to  bright  (grey  level  15). 

DIMENSION  NAME(7) ,10(256) ,NO(256; 

INTEGER  A,B 

ACCEPT"Name  of  image  to  be  examined?  " 
READ(11 ,30) 

30  F0RMAT( SI  3) 

CALL  DUNPACK( NAME, "TEMP") 

CALL  OPEN  (1 ,"TEMP",1 ,IER) 

IF(IER.NE.I) TYPE "OPEN 1  ERRORr  ",IER 
CALL  CFILW  ( "NEW", 2, IER) 

IF(IER.NE.I) TYPE "CREATE  ERR0R=  ",IER 
CALL  OPEN  (2, "NEW", 3 ,IER) 

IF(IER.NE.I) TYPE"OPEN2  ERR0R=  ",IER 
ACCEPT"  Number  of  threshold  value?  ",NUM1 
DO  1  1=0,255 

CALL  RDBLK  ( 1 , I , 10 , 1 , IER ) 

IF(IER.NE.1 )TYPE"RDBLK  ERRORr  ",IER 
DO  2  J  =  1  ,256 
A=IO( J) 

IF( A.LT. NUM1 )G0  TO  3 
NO( J) =15 
GO  TO  2 

3  N0( J) =A 

2  CONTINUE 

CALL  WRBLK  ( 2 , I , NO , 1 , IER) 

IF( IER.NE. 1 )TYPE"WRBLK  ERROR=  ",IER 
1  CONTINUE 

CALL  RESET 

CALL  TREPACK( "NEW") 

ACCEPT"  Name  of  output?  " 

READ( 11 ,30)  NAME(1) 

RENAME  "TPICTN" , NAME 
CALL  DFILW  ("TEMP", IER) 

CALL  DFILW  ("NEW", IER) 

STOP 

END 
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